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ABSTRACT 



A two-level graphical -numerical prediction model is extended so as 
to incorporate some orographical and non-adiabat ic warming or cooling 
effects in a manner such that almost no extra work is added. 

The model is applied to several cases where non-adiabat ic heat 
exchanges are pronounced. The resulting prognoses show a material improve— 
ment over a similar model which omits the non-adiabat ic influence. 
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1, Introduction 



It is a -^eli-kno*n fact that during the winter season cold cP air 
flowing out from continent over the warm ocean surface will acquire heat 
energy from below so that the lower layers of the atmosphere are i#artned. 

As Craddock [l] measured in the case when the arctic air flowed out from 
Iceland to the British Isles, the heating rate may be as large as 65/ly 
per hour and in the course of a day the 1000-500 mb layer thickness could 
increase approximately 800 feet. Hence, in such cases, an appreciable 
error will be produced in the results of those numerical prediction models 
which assume that the air raot^ion is adiabatic. Several techniques have 
been introduced to eliminate this non-adiabat ic error; Reed [5] has 
incorporated Burke's semi-empirical treatment for the air mass transforma- 
tion cP-mP into a graphical model of numerical prognosis. His results show 
considerable improvement over a model omitting the non-adiabat ic heating. 

In the present investigation, the non-adiabat ic heating (or cooling) 
is approximated in a rather crude but extremely simple manner. As a 
result there is no addit ional* work when non-adiabat ic heating and the 
orographic influences [3] are included. 
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2. Development of the prediction equation 



The basic vorticity equation in the (x,y,p,t) coordinates takes the 



form 






( 1 ) 



where is the individual rate-of -change operator following the isobaric 

projection of the three-dimensional motion, 

Jr ttie relative vorticity, 

^ Coriolis force, 

del operator. 

Both J" and refer to values measured on an isobaric surface. 

In Eq, (1) the terms involving the vortex tube and vertical advection 
of vorticity have been neglected, and also the term — ■ V 
Expanding the left side of Eq . (1), 



=:-/\57>\y 



( 2 ) 



By equation of continuity, 

V- V = — 

Eq . (2) becomes 
where 

Assume now that co has a vertical parabolic distribution plus the 
terrain-induced vertical motion. 
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( 4 ) 











are the surface pressure, density and wind 



velocity; H is the height of the terrain; is an upper -pressure level 



constant of proportionality. 

The term on the right side of Eq, (4) involving hj represents the 
terrain induced vertical motion. Here the surface vertical velocity is 
assumed to be proportional to the horizontal wind and the gradient of 



Finally, the terrain-generated vertical velocity is assumed to 
decrease linearly with pressure, and thus we arrive at the first terra constl 
tut ing Ca) in Eq . (4) , 

The second term of Eq . (4) represents the commonly used parabolic dis- 
tribution of vertical velocity assumed to be typical of large-scale press- 
ure systems [4]. Where terrain effects are negligible, this term corre- 
sponds to low-level "isobaric ' velocity divergence and high level con- 
vergence (or vice versa), together with the level of non-divergence at 
Differentiating Eq . (4) with respect to ^ , substituting it into 




terrain height 




• Moreover, 





^ follows that 




Eq . (3) and evaluating the result at 
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( 5 ) 



The next step consists of utilizing the first law of thermodynamics 
to provide a second equation involving CO^y^ • The non-adiabat ic heat 
addition to the surface-500 mb layer will be assumed to be proportional 
to the rate at which the surface air parcels are exposed to a changing 
temperature of the underlying surface. 

Expanding the individual change of the potential temperature, 0 , 

and dividing through by ^ , 



9 dt e 



/ 









( 6 ) 



The subscript p at the ^ operator will be omitted henceforth. 
It can be shown from Poisson's equation 

0 = 

that with pressure held constant 



where is the specific volume and R, the universal gas constant. Sub- 
stituting these identities into Eq • (6) yields 



Bat o( it <y. 



Q 



(7) 



From the hydrostatic equation, 



/v' — _ 

^ a jiy 

substitute into Eq . (7) after multiplying through by ^ , giving 
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It follows that 








( 8 ) 



Equation (8) states that the local change of the rate of change of the 
pressure height with respect to pressure is the result of: a) advection, 

b) vertical motion, and c) the non-adiabat ic heating process represented 
by the last term. 

At this stage, the non-adiabat ic addition of heat will be assumed 
proportional to the expression ^ ' V 7^ where is the wind at the 
surface and is the surface temperature of underlying ocean surface. 
Hence Eq . (8) may then be expressed in the form 



where is the constant of proportionality. Strictly speaking, is not 
a constant but may possibly vary with synoptic situations and also var> 



invest igat ion , 

Now substitute Eq . (4) into Eq . (9) and make geoatrophic approximation in 
the advection term. 




(9) 



slowly with time. However, ^ will be assumed constant for the present 




VH + 





By integrating 




assuming to be a constant, we obtain 
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iy'V%) J.!^, 



+ r^c 



? 



Here h is the thickness between 



/^and , ' 



In the last term of the preceding equation, and are independent 
pressure, hence 






where is a new constant. 
Solving Eq. (10) for Ct). , 






and substituting it into Eq. (5) gives 






i t XMS;^y) — ~f -% ' VH 

* 

' - 4 j] 

\ * V ! 



Rearranging , 






- \/ • <57 1£ ■jxtl’. ^ >y 

(5 /r V/;* 4-1 



n 
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( 10 ) 



( 11 ) 
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Nov, we make the geostrophic approximation for J in Eq . (11). 



1 = 



cri)' 4 



3i r. J 



According to Fjortoft, the horizontal variation of -rr- may be neglected 

•/ / 



JJL 



in comparison to ( . With multiplication by . — — , ve obtain 



lrs^2 1 






( 12 ) 



The Coriolis force term on the right side of Eq • (12) has been 
manipulated as follows: 

rd£ sJ 3,^ - „ 

i 



Now let 



d'ji} 



ir 



/VI • 






r ^ —J-'SI-L ^ 

Therefore Eq. (12) becomes 
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( 13 ) 



■'-Gr-i-'F+Ni:] 



The functionyV'l^varies very slowly with time and may be assumed 
to be a constant for the normal forecast period. The functions G and 
F depend on geographical position only and are invariant with time, 
hence they may be inserted on the left side of Eq. (13). Equation (13) 
shows that the quantity ( ^^—2^ + C^(\ '-t f P* 'f- ) i-S con- 

servative with respect to wind field at /^ . By using the geostrophic 
approximation for^ , i.e. 



A 






^ ^ f 



Eq . (13) may be written in Jacobian form as 



.j 



( 14 ) 



In order to arrive at a more conservative field for adverting the 
conservative property, the Jacobian identity, J(x,y) « J(x+y,y), may be 
applied to Eq . (14), giving 

=- +Gf + (3 + 2„+ ^ + r*/V;r j 05) 

The use of this form implies that the / C-k^ (Sf F J field 

is expected to be more conservative than the field for advection pur- 
poses. This is not necessarily alvays the case since the thickness field 
appears in the advectlng quantity. Nevertheless, in the examples tested, 
Eq . (15) gave good results. 
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To depict the local change of 4- Cf ~t -tWT^J 

shall advect this function for a short cine interval, say twelve hours 

or twenty-four hours. Denote this change by ^ . Since G, P and NT are 

s 

not functions of Cine, 

-t-FtA/Ts] 



~ ^ [^t ] =r (U) 

Moreover, 

•e 

= WC ~(l+C)Z, 1-CZ, 1 =r- A 



Taking Che space average over Eq. (17) and substituting back to Eq. (17), 
yields 






— 



l-HZ 



Lb 






1 






a-K)' 






Averaging n tines. 






7-^ B 



(18) 
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However, this series converges rapidly; and, as shown by FJortoft [2], 
it can be approximated as 

6 S (19) 

In the cases tested in this investigation, it was found that the 
second term on the right side of Eq . (19) can be neglected except when 
a Marked change is indicated. In the routine work where a time saving is 
of great importance, ^ ^ used without much sacrifice in 

accuracy. In utili^ing the high speed electronic computer, Eq . (16) 
could be solved to any desired accuracy. 
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3. Procedure 



In order to provide a comparison to Reed*s results [5], the same 
case has been chosen for prognostication, T^o other series were also 
chosen for testing purposes. One series has the storm located in the 
Gulf of Alaska. Each series contains three consecutive periods with a 
time interval of twelve hours. Except in the same case as Reed's, where 
the prognosis without heating effect was already available, two prog- 
noses for each case have been obtained, one omitting and the other includ- 
ing the heating effect. 

•* 

By omitting the heating effect, Eq • (15) becomes 

~~ -f 01. + ^ i-F ] ( 20 ) 



The following steps are carried out in each prognostication: 

(1). A space-mean chart for is graphically constructed using a 
grid distance of 600 km. 



(2) . With C s 0.5, the current thickness chart is multiplied by 

1/2 (this is achieved by simply tracing every other line on the thickness 

chart) and added to Z . 

o 

(3) . Next the result of step (2) is added to the (G + F) function 
provided by Haltiner and Hesse (3). 

(4) . With N *— 0 / over the ocean and N * -0.05 over the 

land area of Alaska (see Section 4), a (NT ) function chart was constructed 

s 

based on the monthly mean ocean surface Isotherms published by the Hydro- 

graphic Office and the current reported temperatures over the land. 

(Fig. 19). This chart of f(NT ) -'as added to the result of (3). Note 

s 

this step is unnecessary if non-adiabatic heating is omitted. Also it is 
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suggested, since monthly mean surface temperature isotherms over the 
ocean are used, a chart of (G 4 F 4 NT^ ocean) function may be prepared 
for each winter month. This will combine step (4) into step (3) if the 
(NT^) function over land is also approximated by the monthly mean iso- 
therms over the land. Up to this step, a chart of function of 

(Z +Ch 4 'G 4 *F 4 -NT) has been obtained and will be referred to as 
o s 

steering chart henceforth. 

(5) . The chart is subtracted from steering chart to obtain the 
quantity to be advected. 

( 6 ) . Next advect the quantity obtained in step (5) with the steering 
chart for twelve-hour intervals. 

(7) . Subtract the result of step ( 6 ) from the result of step (5) to 
get A. 

( 8 ) . Add 1/2^Z^ to A and multiply the result by 2/3 to obtain B. 

w 

(9) . Average B to get B and multiply by 4/3. 

(10) . Add the result of (9) to ( 8 ), to getAZ . 

o 

( 11 ) . Add A Z to the initial 1000 -mb chart to give the prognostic 

o 

1000 -mb chart of twelve hours later. 
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4. Discussion 



In this investigation, three constants are of major importance, 
namely; C, kj, and N. Ol these three constants, has been discussed in 
the Haltiner and Hesse paper [3], and only C and N will be discussed here. 

(1) . The constant C. 

Several values of C have been used by various investigators. 

In the present investigation, the value of C * 0.5 fits best in the cases 
tested. Values lower than this will give an advection field too weak, 
while higher values give an advecting field too strong. 

(2) . The constant N. ^ 

In this investigation the primary interest is in the area of 
outflow of cP air from inland Alaska to the sea. Over the ocean the value 



of constant N was found empirically to be -0.1 



In the non- 



adiabatic modification of an air mass, the processes which actually take 
place are not simple. The degree and rate of the modification are depend- 
ent upon many factors such as: a) the temperature difference between the 

cP air mass and sea surface; b) the roughness of the sea surface; c) the 
wind speed over the sea surface; d) the length of time the air mass has 
been over the ocean; e) the stability of the air mass; f) the moisture con- 
tent of the air mass; etc. The inter-relations among these factors are not 
well known quantitatively; consequently, it is difficult to determine in a 
specific case the relative importance of the various factors. For example, 
air moving with a higher speed is not necessarily being modified pro- 
portionately more than air moving at a much lower speed. The reason for 
this is that the modification is not only a function of the distance over 
which the air has traveled from the coast, but also a function of the time 
the air mass has been over the ocean. The expression implies 

that when the value of V is doubled, the amount of modification will 

0 
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double, svhich is not necessarily true. Thus, the determination of the 
value of N is a matter of compromise of the many different influences. 

For the relatively small sample tested, the value of N * -0.1 appears to 
be the best choice thus far. However, application of this value shows 
considerable improvement over the other models which treat the non- 
adiabatic heating effect in a different way. 

In the course of this investigation, it appeared that the advecting 
field tends to move the high pressure cell down stream over inland 
Alaska, with the result that the prognostic pressure height was too high 
over the southwestern side of the High and too low over the southeastern 
side of the High. For example, in the case 1500 GMT January 4, 1956 to 
0300 GMT January 5, the high was built up on the east side. In this situa- 
tion, an attempt was made to apply the non-adiabat ic effect over Inland 
Alaska by putting N = -0.05, obtaining a f(NT ) based on the Isotherms on 
the current surface map. Considerable improvement resulted. Such an 
approach provides a continuous temperature field from the inland out to the 
sea. This value of N was also applicable to the other inland cases. The 
combined functions (NTg) for the land and the sea are shown in Figs. 19, 

20, 21, 22, From these figures, it seems feasible to use the monthly mean 
temperature over Alaska land area as the basis of f(NTg) over land, so that 
step (4) can be entirely incorporated into step (3) by preparing a single 
chart of j{0 + F + NTs^ each winter month. Thus a considerable amount 

of time can be saved in daily routine work. Since the non-adiabat ic cool- 
ing process over the land, probably mostly due to radiational cooling, is 
beyond the scope of this investigation, no further discussion will be given. 

In addition to the empirical determination of constants, another 
aspect at the advecting stage in the procedure warrants some discussion. 

As the air mass moves southward over an increasing temperature field, the 
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air mass \vill be N^armed up from belov,/. An assumption was made that the 
air parcel directly in contact with the sea surface will immediately 
acquire the same temperature as that of the sea surface. In this case, 
the air mass tends to be unstable. This instability will enhance the 
upward heat transport, and so increase the modification. On the contrary, 
when the air mass moves northward over a decreasing temperature field, 
by the same assumption it tends to be stabilized. In the first case, the 
thickness of the lower layer will increase and in the latter case the 
thickness will decrease. However the magnitude in the first case of 
increasing thickness is normaTly much larger than the magnitude in the 
latter case of decreasing thickness. Because of this consideration, an 
attempt was made to advect less distance in cold advection and to advect 
the full distance in warm advection. However, this did not improve the 
results, since the property to be advected is a combination of Z^, 

Ch, G, F, and NT^, and not a function of NT^ alone. As a matter of fact, 
since the value of N has been determined empirically, it represents a 
compromise of all the various influences. 
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5. Results 



The results of the tests are displayed in Table I which shows the com- 
parison of the series, with and without non-adiabat ic warming or cooling. 
The statistical parameters computed were based on a seventy -point grid of 
which forty-seven points were over the sea. These grid points are the 
intersections of the meridians and latitude parallels which are 5^ apart 
extending from longitude 165^W to longitude 120^W and from latitude 65^N 
down to latitude 30^N. Due to the convergence of meridians, the points 
are denser over the northern part of the area. Since the area of denser 
points is inland and coastal Alaska, which is a more critical area, the 
values in Table I are more conservative and do not exaggerate the improve- 
ment due to the inclusion of the non-adiabatic heating. 
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Figs. 1 to 16 indicate the results of four cases. These are: 

Case 1) 1500 GMT January 4, 1956 to 0300 GMT January 5, 1956; Cases 2, 

3, 4) 0000 GMT January 13, 1959 to 1200 GMT January 14, 1959, consisting 
of three consecutive prognoses. The other three cases from 0000 GMT 
January 13, 1959 to 1200 GMT January 14, 1959 are not sho'-m. It has also 
been found that the addition of terra NT^ gives no improvement when a 
cyclone comes from the central North Pacific area; however, neither is 
there a decrease in accuracy when f(NTg) is included. 

It can be seen from the figures that the greatest improvement in the 
prognostic pressure heights is found over the area on the west of the storm 
which is the normal location of the cP air. This was true in all the cases 
tested . 

On the basis of the sample tested it may be concluded that this model 
shows very good promise. However, further testing would be desirable to 
determine the best values of the various constants. The non-adiabat ic 
heating terra can readily be incorporated into the programming of baroclinic 
models for the electronic computer on which further testing would be 
greatly facilitated. There is obviously a considerable saving in time 
compared to the Reed method when this model is applied manually in the 
graphical technique . 
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Fi 5 . 1, 1000-mb chart for 1500 GMT 4 January 1956. 
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Fig. 2. 1000-mb prognostic chart for 0300 GMT 5 January 
1956, without heating effect. After Reed* 
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Fig. 3. Actual lOOG-mb chart for 0300 GMT 5 January 1956. 
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Fig. 4. Prog, minus actual. (Fig. 2 minus Fig. 3). 
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lOOO-mb prosnostlc chart for O 3 OO GMT 5 January 

1956, with heating effect. 
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Fig. 6. Prog, with heating effect, minus actual (Fig. 
minus' Fig. 3)* 




t-* 

th 

•H 



25 



1000-mb chart for 0000 GMT 6 January 1959* 
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F10. 8. 1000-mb prosnoatlc chart for 1200 GMT 6 January 
1959 , without heatins effect. 
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1000-mb prognostic chart for 1200 GMT 6 January 
1959 , with heating effect. 
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Fig. 10. Actual 1000-mb chart for 1200 GMT 6 January 1959. 
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1000-mb pro 5 nostlc chart for 0000 GMT 7 January 

1Q59. without heatinff effect, ^ 
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Fig. 12, 1000-mb prognostic chart for 0000 GMT 7 January 
1959 i with heating effect. 
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Fig. 15 . Actual 1000-mb chart for 0000 GHT 7 January 1959. 
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-mb prognostic chart for 1200 GMTTJanuary 
, without heating effect. 
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Fig. 15 . 1000-mb prognostic chart for 1200 GMT January 
1959 * with heating effect. 
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Fig. 16. Actual 1000-mb chart for 1200 GMT 7 January 1959. 
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Fig. 20. Function (N^)ln units of 100 feet Included land area 
prepared at 0000 GMT 6 January 1959 for the prognostication 

of the subsequent twelev hoiu's. 
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. Function (ST.lpwpnred at 1200 GMT 6 January 1959. 
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Fig. 22, Function (ITJ^)prepared at 0000 OMT 7 January 1959. 
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